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摘 要 : 为 探 明 西北 干旱 区 土壤 水 力 参 数 的 空 


间 分 布 特征 ,在 新 疆 地 区 布设 107 个 样 点 ,获取 表层 (0 ~ 5 cm) Ks 


(土壤 饱和 导 水 率 ) FC C 田间 持 水 量 ) 和 PWP( 永 久 著 汞 系数 ) 数 据 ,采用 经 典 统计 学 和 地 统计 学 分 析 土 壤 水 力 参 
数 的 空间 变异 特征 ,并 利用 逐步 回归 和 传递 函数 方程 对 Ks FC 和 PWP 进行 模拟 。 结 果 表 明 :新疆 地 区 Ks FC 和 
PWP 平均 值 分 别 为 10.999 mm - h^! .0.162g.g $80.077 g - o ,变异 系数 为 39.88% ~ 96.07% , 均 表 现 为 中 等 


程度 变异 。 在 区 域 尺度 上 Ks FC 和 PWP 具有 强 


递 函数 方程 相 比 ,多 元 逐步 回归 方程 预测 土壤 水 力 参 数 精度 更 高 , 且 方程 使 用 更 便利 。 


空间 依赖 性 , 变 程 为 97 ~291 km。 与 经 自 变 量 数据 转换 构建 的 传 
基于 土地 利用 、 容 重 、 土 壤 


质地 、 有 机 碳 含 量 和 坡 向 的 多 元 逐步 回归 方程 ,预测 Ks .FC 和 PWP 的 决定 系数 R 分 别 为 0.290、0.494 和 0.491， 


均 方 根 误差 (RMSE ) 分 别 为 2.540 mm . h^! 


,0.039 g + g7' FN 0. 023 g g^, 
土壤 水 力 参 数 的 快速 估算 , VD TY A ZO T A MY EE ATE: AS 7K SCH HE EES XK, 


这 些 方程 的 建立 ,有 利于 西北 干旱 区 


关键 词 : 土壤 饱和 导 水 率 ; 田间 持 水 量 ; KARR RRL; 


土壤 饱和 导 水 率 (soil saturated hydraulic con- 
ductivity, Ks) 、 田 间 持 水 量 (field capacity, FC) 和 永 
ANS IS FB permanent wilting point, PWP) 是 土壤 
重要 的 水 力 参数 。Ks 影响 着 降水 水 分 人 渗 和 产 流 
模式 ,进而 影响 区 域 水 分 再 分 配 过 程 ” 77, FC 是 水 
分 饱和 的 土壤 经 重力 排水 后 土壤 所 能 保持 的 最 大 含 
水 量 , 是 多 数 植物 可 利用 土壤 水 的 上 限 *”-”。PWP 
是 植物 产生 永久 萎 莫 时 的 土壤 含水 量 ,决定 了 植物 
可 利用 土壤 水 的 下 限 ”“”。 同 时 Ks, FC 和 PWP 也 
是 生态 水 文 模型 中 非常 重要 的 敏感 性 参数 ,在 农业 
生产 实践 中 决定 了 农田 需 水 量 与 灌溉 制度 -9 。 相 


关 研究 表明 , Ks, FC 和 PWP 具有 较 强 的 空间 异 质 
性 ”“” 。 因 此 ,准确 揭示 土壤 水 力 参 数 的 空间 分 布 


特征 在 农业 .工程 以 及 水 利 方面 具有 重要 的 意义 。 

目前 ,已 有 学 者 对 土壤 水 力 参 数 的 空间 异 质 性 
进行 了 大 量 研究 。 但 是 由 于 取样 困难 ,测定 时 间 长 
和 耗费 人 力 物 力 较 高 等 原因 ,以 往 土 壤 水 力 参 数 的 
测定 研究 多 集中 于 小 尺度 。 例 如 : 赵 春 雷 等 ”在 黄 


区 域 尺度 ; 


多 元 逐步 回归 ; 新 疆 


坡 面 尺度 表层 Ks WASTE EE” 在 科 尔 
沁 沙 地 对 比 研究 了 不 同 试验 方法 测定 FC 的 差异 。 
Qiao 等 5 在 点 尺度 上 研究 了 黄土 高 原 关 键 带 深 剖 
面 FC 和 PWP 垂直 分 布 特征 。 在 区 域 尺度 上 ,Wang 
等 中 和 Zhao 等 中 采用 网 格 布 点 取样 的 方式 研究 了 
黄土 高 原 区 域 尺度 Ks 的 空间 分 布 特征 。 基 于 遥感 
水 分 数据 , 杨 绍 钵 等 "估算 了 华北 及 东北 地 区 的 土 
HE FC。 以 往 区 域 尺度 土壤 水 力 参 数 的 研究 大 多 集 
中 于 我 国 东部 和 中 部 地 区 ,在 我 国 西北 地 区 则 研究 
较 少 。 

传递 函数 模型 可 以 根据 较 容易 获得 的 土壤 理化 

质 快 速 推算 出 不 易 实地 测量 的 土壤 参数 ,因而 被 
et 5 
Zhao 等 中 基于 土壤 容重 .质地 和 土壤 饱和 含水 量 的 
传递 函数 方程 ,预测 了 Ks 在 整个 黄土 高 原 的 分 布 特 
征 。 利 用 土壤 容重 、 质 地 和 有 机 碳 的 传递 函数 , Dai 
等 5 制作 了 我 国 大 陆地 区 30" x 30” Ks FC 和 PWP 
等 土壤 水 力 参 数 分 布 图 。 传 递 函 数 的 构建 需要 以 大 


土 高 原 北 部 六 道 沟 小 流域 研究 了 坡 面 尺 度 Ks 的 空 
间 异 质 性 。 张 川 等 在 桂 西 北 喀斯特 地 区 研究 了 
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量 实测 的 土壤 水 力 参数 数据 作为 基础 。 然 而 ,有 关 
我 国 西北 干旱 区 土壤 水 力 参 数 的 研究 中 ,存在 研究 
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往 具 有 较 强 的 区 域 适 用 性 ,在 其 他 地 区 建立 的 传递 
函数 很 难 在 西北 干旱 区 获得 较 好 的 预测 结果 ""。 
因此 ,通过 高 密度 的 野外 实地 采样 ,获取 大 量 、 可 靠 
的 实测 数据 ,对 于 认识 西北 干旱 区 土壤 特征 分 布 格 
局 .构建 适宜 本 地 区 的 土壤 水 力 参 数 预 测 方程 是 十 
分 必要 的 。 

本 研究 选取 我 国 新 疆 地 区 ,采用 网 格 布点 的 方 
式 , 实 地 测量 了 新 疆 地 区 土壤 表层 0 ~5 cm Ks、FC 
和 PWP。 采 用 经 典 统计 学 和 地 统计 学 理论 分 析 了 
土壤 水 力 参 数 在 新 疆 地 区 的 空间 分 布 特征 , 且 基 于 
土壤 容重 、 质 地 、 有 机 碳 等 环境 变量 构建 了 预测 Ks. 
FC 和 PWP 的 传递 函数 方程 ,以 期 为 我 国 西 北 干 早 
区 土壤 水 资源 的 科学 管理 .土壤 水 循环 模拟 与 预测 
提供 依据 。 


1 研究 地 区 与 研究 方法 


1.1 研究 区 概况 

研究 区 位 于 我 国 新 疆 (34? 22’ ~ 49?33' N, 
73°22' -96?21'E) (图 1)。“ 三 山 夹 两 倪 ” 构 成 了 新 
性 独特 的 地 形 特点 。 从 南 到 北 分 别 是 昆仑 山 ` 天 山 
和 阿尔 泰山 山脉 ,以 及 中 间 的 塔里木 盆地 和 准噶尔 
例 地 。 中 间 的 天 山 山脉 把 新 疆 地 区 分 为 了 南 疆 与 北 
性 两 部 分 , 南 红 与 北 狂 地 区 分 别 有 塔 克拉 玛 干 沙漠 
与 古 尔 班 通 古 特 沙漠 。 新 疆 地 处 内 陆 远离 海洋 , 气 
候 以 干旱 为 主要 特征 。 北 疆 地 区 多 年 平均 气温 为 
-4 ~9 C ,平均 降水 量 为 150 ~ 200 mm, HERE 
量 为 1 500 ~2 300 mm; 南 疆 地 区 多 年 平均 气温 为 
7 ~14 € ,平均 降水 量 为 25 ~ 100 mm, HERRE 
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Fig.1 Distribution of the sampling sites in Xinjiang 
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为 2 100 ~3 400 mm" 7? , gri Hi DX ^E AS PRE HE 
38 , XUph "^38 , THERE, ERZ UI LUPIS Joc A Se UR 
土 为 主 。 天 然 植被 类 型 主要 有 山地 森林 EE AK 
林山 地 草地 与 荒漠 稀 琉 草原 等 。 
12 样品 采集 与 分 析 

2017 年 8 一 9 月 ,依托 公路 网 ,采用 网 格 法 ( 间 
距 约 为 80 km x80 kmy) 布 点 取样 (图 1)。 选 取 的 采 
样 点 既 考 虑 了 随机 性 ,同时 也 能 代表 采样 点 附近 的 
主要 地 形 特征 和 土地 利用 类 型 。 利 用 GPS 记录 每 
个 点 的 经 纬度 和 海拔 信息 ,使 用 罗盘 仪 测量 样 点 的 
坡度 和 坡 向 。 在 每 个 样 点 用 环 刀 (体积 100 cm ) 采 
集 表层 0 ~5 cm 原状 土壤 样品 1 个 ;在 采样 点 附近 
用 土 钼 采集 表层 0 ~5 cm 扰动 土壤 样品 3 个 并 混合 
为 1 个 扰动 样品 , 土 样 带 回 实验 室 分 析 处 理 。 本 研 
究 共 布设 107 个 样 点 ,由 于 部 分 地 区 土壤 碎 石 含量 
高 ,土壤 原状 土 取样 困难 以 及 没有 道路 到 达 , 因此 ， 
沙漠 和 戈壁 地 区 代表 性 样 点 较 少 。 

采集 的 原状 土壤 样品 用 于 分 析 土 壤 水 力 参 数 。 
Ks 采用 上 供水 恒定 水 头 一 维 土 柱 人 渗 法 路 。FC 与 
PWP 分 别 对 应 土壤 水 势 值 为 -1/3 巴 和 -15 巴 的 
SEA GE 。 土 壤 含 水 量 与 土壤 水 势 之 问 的 关 
系 通过 土壤 水 分 特征 曲线 求 得 ,土壤 水 分 特征 曲线 
采用 离心 机 法 测定 并 用 van Genuchten 模型 拟 合 。 
原状 土 样 经 烘箱 105 % 烘 干 至 恒 重 后 , 称 重 测定 计 
算 土壤 容重 。 扰 动土 壤 样 品 经 风干 后 分 别 过 0. 25 
mm 和 1 mm 网 第 预 处 理 , 用 于 测定 土壤 有 机 碳 含量 
和 土壤 粒 径 组 成 (美国 农业 部 制 ) 。 利 用 重 铬 酸 钾 
氧化 法 测定 土壤 有 机 质 含 量 , 使 用 激光 粒度 仪 测定 
土壤 粒 径 组 成 。 采 样 点 处 的 多 年 平均 降水 量 、 平 均 
气温 和 干燥 度 指 数 等 数据 来 源 于 中 国 科学 院 资源 环 
境 科学 数据 中 心 。 采 样 点 的 土地 利用 组 成 为 : 农 地 
27 个 ,林地 (包括 乔木 和 灌木 林地 )28 个 和 草地 52 
个 (图 1)。 
1.3 ”数据 分 析 方法 

经 典 统计 学 中 主要 采用 变异 系数 (Co，% ) 表示 


空间 变量 的 变异 性 ,其 中 : 


Cv == x100% (1) 
X 


式 中 :5 为 标准 偏差 ;为 平均 值 。 根 据 变异 等 级 分 
Re, Cv<10% 属于 弱 变 异 ,10% < Cv<100% JRF HF 
等 变异 ,Co > 100% JE Fane RY 。 

地 统计 学 空间 变异 理论 主要 基于 半 方 差 函 数 。 


李 祥 东 等 :西北 干旱 区 土壤 水 力 参 数 空间 变异 与 模拟 


ENA yA) 是 计算 变量 在 空 
下 的 样本 方差 值 ,计算 公式 为 : 


空间 位 置 相 隔 一 定 距 离 


y(h) = Sy lel) 2% +h)) (2) 


USE. 
式 中 六 为 滞后 距离 ia 为 样 点 总 数 ;x; 表示 样 点 位 
置 ; Z(%i) 为 随机 变量 Z 在 点 x =xi 处 的 值 ; ZG; + 
有 ) 为 相距 点 x, Jy h 和 远 处 的 值 ;N(h) 为 滞后 距离 是 有 


时 的 样本 对 数 。 半 方差 函数 反映 了 变量 在 空间 上 的 
相关 性 。 通 过 半 方 差 图 可 以 计算 出 空间 变量 的 块 金 


ECC) EKE CC + C) RIZR E (Range) 7, 48 
据 Cambardella 的 划分 标 NEUSS 空间 异 质 比 Cy/ 
(Co+C) 友 25% 时 表示 强 的 空间 依赖 性 ,25% < Cy/ 
(Co +C) <75% 时 表示 中 等 空间 依赖 性 ,CoA( Co + 
C) 二 75% 时 表示 弱 的 空间 依赖 性 ”。 

选择 海拔 Mu) IE VER A BRL AL. 
有 机 碳 含量 .多 年 平均 降雨 量 E35 HCl TERRE TR 
数 和 土地 利用 类 型 作为 环境 变量 ,采用 多 元 逐步 回 
归 方 程 预测 土壤 水 力 参 数 , 其 中 逐步 回归 变量 的 选 
择 与 移 除 标 准 是 依据 变量 引入 后 显著 性 PP 值 0.05 
(引入 ) 与 0.1( 移 除 )。 同 时 选取 与 土壤 水 力 参 数 
显著 相关 的 环境 变量 , 自 变 量 经 对 数 .倒数 SET JT 
方 4 种 常见 的 数据 转换 后 进行 逐步 回归 分 析 , 构 建 
土壤 水 力 参 数 传递 函数 方程 。 土 壤 水 力 参 数 数据 集 
随机 分 为 75% 的 训练 样本 (80 4+) FI 2596 的 验证 样 
本 (27 个 ) 。 利 用 训练 样本 构建 多 元 回归 方程 与 传 
递 函 数 ,利用 验证 样本 对 所 构建 的 土壤 水 力 参 数 预 
测 方 程 进行 精度 检验 。 本 研究 采用 决定 系数 (R ) 
和 均 方 根 误差 (RMSE ) 来 评价 逐步 回归 方程 和 传递 
函数 模型 的 精度 。 
1.4 数据 处 理 

使 用 SPSS 20. 0 进行 各 变量 的 基本 统计 分 析 、 
Pearson 相关 分 析 和 逐步 回归 分 析 , 使 用 GS+9.0 进 
行 土壤 水 力 参 数 半 方 差 函 数 计算 与 拟 合 ,使 用 Arc- 
GIS 10. 6 进行 土壤 水 力 参数 克 里 格 插值 及 制图 。 


2 结果 与 讨论 


2.1 土壤 水 力 参 数 的 基本 统计 特征 

表 1 显示 了 土壤 水 力 参数 的 统计 特征 值 。 新 疆 
地 区 Ks, FC 和 PWP 平均 值 分 别 为 10. 999 mm * 
h~! 0.162 gg-! 和 0.077 g.g-!。 相 比 于 我 国 东 
南部 地 区 ,新疆 地 区 具有 较 高 的 Ks, 这 与 本 地 区 较 


R1 土壤 水 力 参数 撒 述 性 统计 特征 
Tab.1 Descriptive statistic characteristics of 


soil hydraulic parameters 


土壤 水 饱和 导 水 率 间 持 水 量 ES Eo 

力 参 数 /(mm*h^) gtg)  /(g+g7!) 
最 小 值 0.060 0. 048 0.013 
最 大 值 62. 968 0.342 0.174 
平均 值 10.999 0. 162 0.077 
标准 差 10.567 0. 065 0.037 
变异 系数 /9% 96.072 39. 875 48.139 
偏 度 1.873 0. 824 1.007 
峰 度 5.370 0. 604 0.645 
K-S 0.019 0. 548 0. 107 
分 布 类 型 对 数 正 态 正 态 正 态 


高 的 砂粒 含量 和 较 低 的 黏 粒 含量 有 关 ' 玉 ; 。 黄 土 高 
原 地 区 "和 喀斯特 地 区 ”表层 土壤 Ks 分 别 为 
13.90 ~42.00 mm * h^! $/549.17 mm * h^ ! ,高 于 本 
研究 区 域 。 这 可 能 是 因为 本 研究 区 中 沙漠 地 带 无 法 
到 达 , 采 样 点 稀少 ,因而 低估 了 新 疆 地 区 平均 Ks 
其 次 陕 北 黄土 高 原 地 区 由 于 本 身 砂粒 含量 较 高 ”， 
因而 表层 Ks 大 于 新 疆 非 沙 漠 地 区 。Dai 0° 4 
出 ,FC 和 PWP 较 高 值 分 布 在 我 国 东 南部 地 区 ,西北 
干旱 半 干 旱地 区 FC 和 PWP BEE Wu 等 (中 指出 ， 
我 国 最 低 的 FC 在 新 疆 南 部 塔克拉玛干 沙漠 地 区 与 
甘肃 北部 阿拉 善 沙漠 地 区 ,西北 和 干旱 区 较 低 的 FC 
与 本 地 区 较 高 的 砂粒 含量 有 关 。FC 和 PWP 在 美国 
干旱 地 区 中 分 别 为 0.210 g+ g' Fl0. 100 g- g^, 
在 黄土 高 原 地 区 "分 别 为 0.228 g- g^ 和 0. 130 
T , 均 显 著 高 于 本 研究 地 区 ,这 也 与 西北 地 区 较 


g'g 
高 的 砂粒 含量 有 关 。 土 壤 水 力 参 数 Ks FC 和 PWP 


的 Cv Æ 10% ~90% ,为 中 等 程度 变异 。 区 域 尺度 
上 黄土 高 原 地 区 土壤 表层 Ks 的 Cv 为 206% ,为 强 
空间 变异 ") ,FC 和 PWP 为 中 等 程度 变异 ”" 。 通 
过 偏 度 、 峰 度 以 及 单 样本 K-S 检验 结果 表明 ,FC 和 
PWP 为 正 态 分 布 , 而 Ks 经 底数 为 10 的 对 数 变 换 后 
满足 正 态 分 布 。 
2.2 ”基于 地 统计 学 的 土壤 水 力 参 数 空间 变异 特征 
及 分 布 格局 

利用 地 统计 学 半 方 差 理 论 对 土壤 水 力 参 数 的 空 
间 变 异性 进行 分 析 ,结果 见 表 2。Ks 与 PWP 的 最 佳 
拟 合 模型 为 高 斯 模型 , FC. 最 佳 拟 合 模型 为 指数 模 
型 。 土 壤 水 力 参 数 具 有 较 小 的 块 金 值 (0. 001 ~ 
0.051) ,表明 随机 因素 引起 的 土壤 水 力 参 数 变异 性 
较 小 。Ks FC 和 PWP 空间 异 质 比分 别 为 16. 90% 、 
0.20% 和 23.30% , 均 小 于 25% ,表明 在 新 疆 地 区 土 
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R2 土壤 水 力 参数 半 方 差 函 数 模 型 及 其 结构 参数 


Tab.2 Semivariance model and structural parameters of soil hydraulic parameters 


土壤 水 力 参 数 模型 块 金 值 基 台 值 空间 异 质 比 /% 变 程 /km 决定 系数 残 差 平方 和 
饱和 导 水 率 高 斯 模型 0.051 0. 305 16.90 97.00 0.555 0.014 
间 持 水 量 * 指数 模型 0.001 0.431 0.20 291.00 0. 788 2.770 
BGA 高 斯 模型 0.031 0. 135 23.30 138. 60 0.753 0. 182 


注 :a 表示 数据 扩大 了 10 fits 


壤 水 力 参 数 具 有 强 的 空间 依赖 性 。 变 程 反映 空间 变 
量 在 多 远 距 离 范 围 内 具有 空间 相关 性 。Ks $697. 00 
km 范围 内 具有 空间 依赖 性 。FC 与 PWP 变 程 分 别 
为 291.00 km 和 138.60 km, Ks 在 黄土 高 原 区 域 尺 
度 空间 异 质 比 为 48. 6096 为 中 等 空间 依赖 性 , 变 程 
24 68 km? 。 表 明 土 壤 水 力 参 数 在 新 疆 地 区 具有 更 
大 的 空间 相关 性 。 邵 明 安 等 中 指出 ,最 佳 采样 点 间 
距 应 小 于 变 程 的 一 半 , 因 此 ,Ks 5j PWP 应 缩小 采样 
间距 ,增加 采样 点 ,而 FC 可 以 适当 增加 采样 间距 ， 
减少 采样 点 数据 。 

半 方 差分 析 表 明 ,Ks、FC 和 PWP 在 新 疆 地 区 具 
有 明显 的 空间 结构 ,为 直观 反映 土壤 水 力 参 数 在 新 
己 地 区 的 空间 分 布 格局 ,利用 基于 半 方 差 理 论 的 普 
通 克 里 格 插 值 法 对 土壤 水 力 参 数 进行 插值 作 图 (图 
2)。 因 新 疆 沙漠 地 区 以 及 东部 戈壁 地 区 采样 点 较 
少 ,插值 结果 不 具有 代表 性 ,因此 ,对 这 些 地 区 不 考 
虑 。 天 山东 部 地 区 .昆仑 山西 部 地 区 和 阿尔 泰山 北 
部 地 区 具有 较 高 的 Ks ;塔里木 盆地 边缘 ET AR d 
地 边缘 和 天 山中 部 地 区 具有 较 低 的 Ks。 新 疆 西 部 
和 塔里木 河流 域 中 上 游 土 壤 具 有 较 高 的 FC 与 
PWP; 伊 犁 河谷 地 区 以 及 新 疆 东 北部 区 域 具 有 较 低 
的 FC 和 PWP 值 。 本 研究 中 新 疆 地 区 土壤 水 力 参 
数 分 布 格局 与 Dai 等 "研究 相 比 ,在 高 值 分 布 区 相 


(2) 饱和 导 水 率 /mm + hi?) 


B <0.10 


0.13-0.16 >02 
图 2 新 疆 地 区 土壤 水 力 参 数 空间 分 布 


Fig.2 Spatial distribution of soil hydraulic parameters in Xinjiang 


(b) 田 间 持 水 量 /(g g» 
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对 一 致 ,而 在 低 值 分 布 区 略 有 差异 。 土 壤 水 力 特征 
高 值 分 布 区 域 往往 伴 有 较 高 的 植被 盖 度 ,因此 ,研究 
土壤 水 力 特 征 有 利于 加 深 西北 干旱 区 植被 水 分 利用 
过 程 的 认识 ™。 
2.3 土壤 水 力 参数 与 环境 变量 的 相关 关系 

为 探究 环境 因子 对 土壤 水 力 参 数 的 影响 ,利用 
Pearson 相关 系数 分 析 环 境 变 量 与 土壤 水 力 参 数 的 
相关 性 ( 表 3) Pearson 相关 分 析 表 明 : Ks 与 海拔 、 
坡 向 坡度 和 砂粒 含量 呈 显 著 正 相关 (P<0.05 ) ,与 
土壤 容重 AM , 粉 粒 含量 和 多 年 平均 气温 呈 显 著 负 
相关 (P<0.05)。 人 砂粒 含量 是 影响 土壤 饱和 导 水 的 
最 显著 因子 ,Pearson 相关 系数 为 0.426。Ks 与 土壤 
有 机 碳 含量 .多 年 平均 降水 和 和 干燥 度 指数 没有 明显 
的 相关 关系 (P >0.05)。FC 和 PWP 与 环境 因子 具 
有 相似 的 相关 关系 。 坡 向 坡度、 容重 .土壤 质地 
(BER AS RL A DE) 和 有 机 碳 含量 显 著 影 响 土 壤 
FC 和 PWP(P <0.01)。 土壤 容重 与 FC H PWP 呈 
显著 负 相 关 , Pearson 相关 系数 最 大 ,为 -0.686 和 
-0.556。 和 气象 因子 对 FC 与 PWP 没有 明显 的 相关 
性 (已 >0.05 ) 。 这 与 前 人 研究 相 一 致 ,土壤 容重 和 
质地 对 土壤 水 力 参 数 具 有 显著 相关 关系 和 7) Ks 
与 土地 利用 类 型 显著 相关 (P <0.01), 单 因素 方差 
分 析 (ANOVA) 表 明 , 农 地 Ks 平均 值 为 5.313 mm - 


(28 € SK (e + g 


EN-004 [1008-010 
EE004-0.06 [E89 0.10-0.12 
[006-008 BN > 0.12 


李 祥 东 等 :西北 干旱 区 土壤 水 力 参 数 空 


间 变 异 与 模拟 


R3 环境 变量 描述 及 其 与 土壤 水 力 参 数 的 Pearson 相关 关系 
Tab.3 Descriptions of environmental variables and the Pearson’s correlation coefficients between soil hydraulic parameters 
变量 平均 值 标准 差 变异 系数 /% SS 
饱和 导 水 率 田间 持 水 量 BEAM 
海拔 * 1 171.92 627.22 53.52 0.217* 0.115 0. 084 
坡 向 22.67 66.50 293. 32 0.305 * 0.198 六 0.240 * 
坡度 1.68 5.37 319.01 0.239 * 0. 306 ** 0.351 
容重 1.33 0.24 17.73 -0.333 * -0.686 * -0.556 * 
Ae t 16.17 8.45 52.26 -0.382 *™ 0.293 '* 0.279 
粉 粒 含量 45.22 17.78 39.31 -0.320 * 0.587 * 0.575" 
fb yar Se f 38.61 23.45 60. 75 0.426 ** -0.554 * -0.541 * 
有 机 碳 含量 * 9.35 16.92 180. 95 0.077 0.439 ** 0.491 
年 平均 降水 量 142.45 97.87 68.71 0. 158 0. 049 0.155 
多 年 平均 气温 6.30 4.17 66.21 -0.232* -0.020 -0.145 
干燥 度 指 数 " 9.01 6.28 69. 72 -0. 104 0.023 - 0.081 
土地 利用 类 型 i i = 0.441 ** -0.072 0.006 


h ,显著 小 于 林地 和 草地 (已 <0.05)。FC 和 PWP 
与 土地 利用 类 型 不 相关 (P >0.05)。 环境 因子 对 土 
壤 水 力 参 数 的 影响 因 研 究 地 区 和 研究 尺度 的 不 同 而 


注 :a 表示 数据 进行 了 以 10 为 底 的 对 数 转换 ; * 表示 在 0. 05 水 平 ( 双 侧 ) 上 显著 相关 ; 


** 表 示 在 0.01 水 平 ( 双 侧 ) 上 显著 相关 。 


2.4 土壤 水 力 参 数 多 元 
拟 合 
土壤 水 力 参 数 传 统 测 定 方法 费时 费力 ,限制 了 


逐步 回归 与 传递 函数 方程 


有 不 同 的 表现 。 例 如 陕西 六 道 沟 小 流域 坡 面 尺度 ， 
Ks 与 土壤 有 机 碳 含量 没有 显著 相关 关系 (P > 
0.05) ,而 在 黄土 高 原 区 域 尺度 上 则 与 有 机 碳 含 
量 旦 显著 负 相关 ,相关 系数 为 -0.15(P<0.01)"。 
因此 ,人 研究 土壤 水 力 参 数 的 变异 性 ,必须 考虑 地 理 位 


置 与 研究 尺度 特征 。 


表 4 土壤 水 力 参 数 模 拟 的 多 元 线性 逐步 回归 方程 、 传 弟 


区 域 尺度 上 土壤 水 力 参 数 的 相关 研究 。 本 研究 随机 
选取 75% 的 样 点 构建 多 元 逐步 回归 方程 与 传递 函 
数 , 利 用 剩余 25% 的 样 点 数据 对 构建 的 预测 方程 进 
行 精度 验证 ( 表 4 和 图 3)。 构 建 的 逐步 回归 方程 对 
Ks .FC 和 PWP 模拟 的 决定 系数 R 分 别 为 0. 454, 
0.710#10.571,RMSE 4} 3I] 292. 586mm * h^! , 


数 方程 和 模型 预测 精度 


Tab.4 Multiple linear stepwise regression equation and pedotransfer function equation simulated with soil hydraulic 


parameters as well as their prediction accuracies 


训练 组 验证 组 
模拟 方法 土壤 水 力 参 数 方程 pg RMSE /(mm * h^! , ge RMSE /(mm -h E 
gg!) gg) 
逐步 线性 logKs =1.310 +0. 164 x LU -0.964xBD+0.010 x Sand + 0.454" 2.586 0.290 * 2.540 
回归 方程 1.623 x10 ? x SA 
FC 20.408 -0. 165 x BD -0.996 x 10 ^? x Sand + 1.92 x 0.710 * 0.037 0.494 * 0. 039 
10-3 x SOC 
PWP =0. 107 +0. 787 x 10 ^? x Silt -0. 056 x BD + 1.316 x 0.571 * 0.025 0.491 * 0. 023 
107? x SOC 
传递 函数 方程 logKs- -5.854 41.605 x logSand - 3.066 x BD? 40.174 x 0.607 ** 2.251 0.269 ** 2.713 
LU +6. 678 x BD 40.357 x logSOC +0. 023 x Clay 
FC 20.459 —0. 068 x BD? 20.018 x Sand^5 40.017 x 0.740 * 0. 035 0.481 * 0. 039 
SOC? —0.066 x logSilt 
PWP =0. 183 46.231 x10~° x Sil? —0. 134 x BD®5 + 0.621 * 0. 024 0.475 * 0. 024 


注 :LU 为 土地 利用 方式 ;BD 为 容重 ;Sand 为 砂粒 含 


0.010 x SOC?? +0. 182 x Sand ^! 


量 ;SA 为 坡 向 ;SOC 为 有 机 碳 ;Sil S SUBE 


含量 ; ** 表 示 在 0.01 水 于 


FF 上 模型 模拟 显著 。 
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Al3 基于 多 元 线性 


实测 FC/(g * g» 
回归 方程 和 传递 函数 方程 的 验证 组 土壤 水 力 参 数 实测 值 与 模拟 值 对 比 


实测 PWP/(g + g5 


Fig.3 Measured and simulated soil hydraulic parameters based on the datasets validated with multiple linear stepwise regression 


and pedotransfer function equation 


0.037 g.g '#l0.025 g.g 。 构 建 的 传递 函数 方 
程 对 Ks, FC 和 PWP 模拟 的 决定 系数 RO 分 别 为 
0.607 .0.740 和 0. 621, RMSE 分 别 为 2. 251 mm - 
h',0.085g-g 'fi0o.024g* g ', 5j FC fi] PWP 
相 比 ,土壤 Ks 的 预测 精度 较 低 , 说 明 土壤 Ks 受 随 机 
因素 的 影响 更 大 ,更 具有 区 域 独特 性 ,这 与 半 方 差 研 
究 中 Ks 变 程 较 小 的 结果 一 致 ( 表 2)。 与 训练 组 数 
据 相 比 ,验证 组 的 RC 具有 不 同 程度 的 降低 。 与 以 往 


研究 不 同 ”) ,通过 验证 组 数据 发 现 ,逐步 回归 方程 
不 仅 结 构 简 单 且 对 土壤 水 力 参 数 的 预测 精度 更 高 ， 
优 于 自 变 量 经 数据 转化 的 传递 函数 ( 表 4 和 图 3)。 
预测 Ks FC 和 PWP 的 多 元 逐步 线性 回归 方程 验证 
组 决定 系数 RR 分 别 为 0. 290,0. 494 和 0. 491, 
RMSE 分 别 为 2. 540 mm + h^',0.039 g - g 和 
0.023 g' g … 。 土 壤 水 力 参 数 实测 值 与 模拟 值 对 比 
研究 表明 (图 3) ,土壤 Ks 在 较 低 值 时 模拟 精度 较 
高 ,在 高 于 10 mm -h~ 时 模拟 精度 显著 降低 ,在 高 
于 20 mm : h 时 预测 方程 会 低估 Ks。 多 元 逐步 回 
归 与 传递 函数 对 FC 和 PWP 的 模拟 较 一 致 ,预测 方 
EX FC 720.1 ~0.2 g- g` 区 间 以 及 PWP 在 
0.05 ~ 0.10 g - g ”区间 预测 精度 较 高 ,但 在 较 低 值 
时 会 高 估 FC 和 PWP, 在 较 高 值 时 会 低估 FC 和 
PWP, Zhao 等 中 基于 黄土 高 原 数 据 构建 的 传递 方 
程 能 够 解释 Ks 变异 的 31% ,高 于 本 研究 预测 精度 。 
Ks 预测 精度 较 低 与 土壤 Ks 本 身 较 大 的 变异 性 有 
关 , 例 如 在 喀斯特 地 区 "Ks 变异 系数 为 95% ,黄土 
高 原 北 部 坡 面 尺度 变异 系数 为 58% ,黄土 高 原 区 
域 尺 度 中 变异 系数 为 67% ,本 研究 Ks 变异 系数 为 
96.07% , 远 高 于 以 上 地 区 。 此 外 ,考虑 到 模型 的 易 


用 性 ,本 研究 考虑 了 12 个 相对 容易 获取 的 潜在 影响 
Ks 的 因子 。 然 而 ,由 于 Ks 的 影响 因素 众多 且 关 系 
复杂 ,本 研究 所 考虑 的 因素 不 足以 完整 地 揭示 Ks 的 
全 部 变异 特征 ,最终 导致 Ks 的 精度 偏 低 。 王 子 龙 
等 中 研究 表明 ,团聚 体 也 对 Ks 有 影响 ,然而 ,团聚 
体 测定 比较 复杂 ,因此 本 研究 并 未 考虑 。Qiao 等 
a FC 和 PWP 的 估计 精 
度 , 预 测 R 分 别 低 于 0.333 和 0.321 eae 
所 建立 的 预测 方程 。 土 地 利用 类 型 .容重 砂粒 含 
和 坡 向 是 影响 Ks 最 重要 的 因子 。 容 重 、 né 
成 (砂粒 和 粉 粒 ) 和 有 机 碳 含量 是 影响 FC 和 PWP 
最 重要 的 因子 。 


3 结论 


(1) 新 疆 地 区 土壤 水 力 参 数 具 有 明显 的 空 


间 变 


异性 ,变异 程度 为 中 等 。Ks FC 和 PWP 在 本 地 区 具 
有 强 的 空间 依赖 性 , 变 程 为 97 ~291 km, 


(2) 土地 利用 类 型 .容重 和 土壤 质地 是 影响 Ks 
变异 的 最 主要 因子 。 土 壤 容 重 、 质 地 和 有 机 碳 含 量 
对 FC 和 PWP 具有 重要 影响 。 

(3) 与 经 典 多 元 逐步 回归 方程 相 比 ,基于 自 变 
量 经 数据 转化 构建 的 传递 函数 方程 较 复 杂 且 没有 提 
高 模型 的 模拟 精度 。 多 元 逐步 回归 方程 能 够 解释 
Ks 变异 的 29.00% 。 逐 步 回 归 方 程 对 FC 与 PWP 模 
拟 效果 较 好 ,能 够 预测 FC 与 PWP 变异 的 49. 40% 
与 49.10% 。 


参考 文献 ( References : 


(1) 邵 明 安 , 王 全 九 , 黄 明 斌 . 


土壤 物理 学 [ MJ. 北京 :高 等 教育 出 


201911.00068v1 


chinaXiv 


(2) 


(3) 


(4) 


(5] 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


ChinaXiv 合 作 期 刊 


李 祥 东 等 :西北 干旱 区 土壤 水 力 参 数 空间 变异 与 模拟 


版 社 ,2006: 228 - 261. (Shao Ming'an, Wang Quanjiu, Huang 
Mingbin. Soil Physics( M). Beijing: Higher Education Press ,2006 ; 
228 -261. ] 

赵 春 雷 , 邵 明 安 , 贾 小 旭 . 黄土 高 原 北部 坡 面 尺度 土壤 饱和 导 
水 率 分 布 与 模拟 [J]. 水 科学 进展 ,2014,25 (6):806 - 815. 


[Zhao Chunlei ,Shao Ming'an ,Jia Xiaoxu. Distribution and simula- 


tion of saturated soil hydraulic conductivity at a slope of Northern 
Loess Plateau( J]. Advances in Water Science,2014 ,25 (6) :806 
-815. ] 

SORS. 土壤 学 [ Mj. 北京 :中国 农 业 出 版 社 ,2000:60. ( Huang 
Changyong. Pedology( M). Beijing: China Agriculture Press ,2000 ; 
60.) 

Yan N N, Wu B F, Huang H P, et al. Mapping of field capacity and 


wilting coefficient in China for drought monitoring ( C ) // Proceed- 
ings. 2005 IEEE International Geoscience and Remote Sensing 
Symposium. America; IEEE , 2005 , 4 433 — 4 435, doi: 10. 1109/ 
IGARSS. 2005. 1525903. 

Reynolds C A, Jackson T J, Rawls W J. Estimating soil water-hold- 
ing capacities by linking the food and agriculture organization soil 
map of the world with global pedon databases and continuous pedo- 
transfer functions( J]. Water Resources Research ,2000 ,36 (12) ; 
3 653 -3 662. 

Mohanty M, Sinha N K, Painuli D K, et al. Modelling soil water 
contents at field capacity and permanent wilting point using artifi- 
cial neural network for Indian soils( J). National Academy Science 
Letters ,2015 ,38 (5) :373 - 377. 

Wang Y Q,Shao M A,Liu Z P,et al. Regional-scale variation and 
distribution patterns of soil saturated hydraulic conductivities in 
surface and subsurface layers in the loessial soils of China( J). 
Journal of Hydrology ,2013 ,487 :13 —23. 

Zhao C L,Shao M A, Jia X X, et al. Using pedotransfer functions to 
estimate soil hydraulic conductivity in the Loess Plateau of China 
(J). Catena ,2016 ,143 :1 - 6. 

张 川 , 陈 洪 松 , 张 伟 ,等 . 喀斯特 坡 面 表层 土壤 含水 量 、 容 重 和 
饱和 导 水 率 的 空间 变异 特征 [机 . 应 用 生态 学 报 ,2014 ,25 (6): 
1 585 -1 591. ( Zhang Chuan, Chen Hongsong, Zhang Wei, et al. 


Spatial variation characteristics of surface soil water content , bulk 
density and satureated hydraulic conductivity on Karst slope( J). 
Chinese Journal of Applied Ecology, 2014, 25 (6) :1 585 - 
1 591.) 

SHA XER, Bem E , S. 科尔沁 沙丘 - 草 甸 田间 持 水 量 的 
Zi ix 5 SE CJ). 干旱 区 研究 ,2016 ,33(6) :1 336-1 344. 


(Gao Xiaoyan, Liu Tingxi, Duan Limin, et al. Determination and 


effect of dune-meadow field capacity in the Horqin Sandy Land 
(JJ. Arid Zone Research ,2016,33 (6) :1 336 -1 344. J 

Qiao J B,Zhu Y J,Jia X X, et al. Pedotransfer functions for estima- 
ting the field capacity and permanent wilting point in the critical 
zone of the Loess Plateau , China( JJ. Journal of Soils & Sediments, 
2019 ,19(1) :140 - 147. 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 


杨 绍 锣 , 吴 炳 方 , 关 娜 娜 . 基于 AMSR-E 数据 估 测 华北 平原 及 东 
北 地 区 土壤 田间 持 水 量 [J]. 土壤 通报 ,2012,43(2) :301 - 305. 
[Yang Shao’e, Wu Bingfang,Yan Nana. Estimation of soil field ca- 


pacity in north China plain and Northeast China based on AMSR-E 
data( J]. Chinese Journal of Soil Science, 2012,43 (2) :301 - 
305. J 

Dai Y J, Shangguan W , Duan Q Y,et al. Development of a China 
dataset of soil hydraulic parameters using pedotransfer functions for 
land surface modeling[ JJ. Journal of Hydrometeorology , 2015 , 14 
(3) :869 - 887. 

Zheng Z Y, Ma Z G,Li M X,et al. Regional water budgets and 
hydroclimatic trend variations in Xinjiang from 1951 to 2000( JJ. 
Climatic Change,2017 ,144(3) :447 - 460. 

Zhang R P, Liang T C, Guo J, et al. Grassland dynamics in re- 
sponse to climate change and human activities in. Xinjiang from 
2000 to 2014 ( J]. Scientific Reports, 2018, 8:2 888, doi; 10. 
1038/s41598 -018 - 21089 -3. 

Cambardella C A, Moorman T B , Novak J M ,et al. Field-scale vari- 
ability of soil properties in central Iowa soils( JJ. Soilence Society 
of America Journal ,1994 ,58(5) :1 501 - 1 511. 

HÆ AEE, FA, T. 黄土 区 坡 沟 系 统 容重 ,饱和 导 水 率 
和 土壤 含水 量变 化 分 析 [J]. 干旱 区 研究 ,2018 ,35 (2) :315 - 
324. ( Gan Miao, Jia Yuhua, Li Tongchuan ,et al. Variation of soil 


bulk density, saturated hydraulic conductivity and soil moisture 
content in a slope-gully unit on the Northern Loess Plateau[ J]. Ar- 
id Zone Research ,2018 ,35(2) :315 —324. ) 

Wu X T,Lu G H, Wu Z Y ,et al. An integration approach for map- 
ping field capacity of China based on multi-source soil datasets 
(JJ. Water,2018 ,10 :728 , doi: 10. 3390/w10060728. 

Santra P, Kumar M,Kumawat R N,et al. Pedotransfer functions to 
estimate soil water content at field capacity and permanent wilting 
point in hot arid Western India( J). Journal of Earth System Sci- 
ence,2018,127(3) :35 - 50. 

陈 亚 宁 , 李 卫 红 , 陈 亚 鹏 ,等 . 荒漠 河岸 林 建 群 植物 的 水 分 利用 
whee or). 干旱 区 研究 ,2018 ,35 (1) :130 - 136. ( Chen Ya- 


ning,Li Weihong, Chen Yapeng, et al. Water use process of con- 
structive plants in desert riparian forest( J]. Arid Zone Research, 
2018 ,35(1) :130 - 136. J 

于 冬 雪 ,机 小 旭 , 黄 来 明 , 等 . 黄土 区 不 同 土 层 土壤 容重 空间 变 
异 与 模拟 [中 .土壤 学 报 ,2019 ,56(1):55 -64. (Yu Dongxue, 


Jia Xiaoxu , Huang Laiming , et al. Spatial variation of soil bulk den- 


sity in different soil layers in the loess area and simulation( J]. Ac- 
ta Pedologica Sinica ,2019 ,56(1) :55 -64. J 

王子 龙 ,赵勇 钢 , 赵 世 伟 ,等 . 退耕 典型 草地 土壤 饱和 导 水 率 及 
其 影响 因素 研究 [J 了. 草地 学 报 ,2016 ,24(6):1 254 -1 262. 
(Wang Zilong, Zhao Yonggang, Zhao Shiwei, et al. Study on soil 


saturated hydraulic conductivity and its influencing factors in typi- 
cal grassland of farmland conversion ( JJ. Acta Agrestia Sinica, 


2016 ,24(6) :1 254 -1 262. ) 


Ji 
um 
x 
Ex 
E 


Spatial Variability and Simulation of Soil Hydraulic Parameters 
in Arid Northwest China 


LI Xiang-dong'*, SHAO Ming-an'^^, ZHAO Chun-lei? 
(1. State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau „Institute of Soil and Water Conservation , 
Chinese Academy of Sciences and. Ministry of Water Resources , Yangling 712100 , Shaanxi , China ; 
2. College of Natural Resources and. Environment , University of Chinese Academy of Sciences Beijing 100190 , China; 
3. Key Laboratory of Ecosystem Network Observation and Modeling , Institute of Geographic Sciences and 
Natural Resources Research , Chinese Academy of Sciences , Beijing 100101 , China) 


Abstract; Soil saturated hydraulic conductivity ( Ks) , field capacity (FC) and permanent witling percentage 
(PWP) are the important soil hydraulic parameters. These parameters are also the critical factors in bio-hydrologic 
models. In order to understand the regional-scale spatial variability of these hydraulic parameters, the Ks, FC and 
PWP data (0 —5 cm soil layer in depth) were obtained from 107 sampling sites in Xinjiang. Both traditional statis- 
tics and geostatistics were used to explore the spatial variations of Ks, FC and PWP. Estimations of these parameters 
were further conducted by multiple stepwise regression (MSR) and pedotransfer functions ( PTFs) . The results indi- 
cated that the mean values of Ks, FC and PWP were 10.999 mm - h^! ,0.162 g + g | and 0.077 g * g ! respec- 
tively. These three hydraulic parameters varied moderately with Cvs of 39.88% — 96. 0796 , and had strong spatial 
correlations in a range of 97 —291 km in Xinjiang. MSR performed better than PTFs in estimating soil hydraulic pa- 
rameters. Based on the land use , bulk density ,soil texture , soil organic carbon content and slope aspect ,the determi- 
nation coefficients of the MSR models for Ks, FC and PWP were 0. 290 ,0. 494 and 0. 491 , and the root mean square 
errors were 2.540 mm * h^! ,0.039 g - g ' and 0.023 g - g ^! ,respectively. The development of these equations is 
beneficial to the rapid estimation of soil hydraulic parameters and can thus provide the key parameters for the agri- 
cultural irrigation and eco-hydrological models in arid Northwest China. 

Key words: soil saturated hydraulic conductivity; field capacity; permanent witling percentage; regional scale; 


multiple stepwise regression; Xinjiang 


